Anomalous scaling of conductivity in integrable fermion systems 
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We analyze the high- temperature conductivity in one-dimensional integrable models of interacting fermions: 
the t-V model (anisotropic Heisenberg spin chain) and the Hubbard model, at half-filling in the regime cor- 
responding to insulating ground state. A microcanonical Lanczos method study for finite size systems reveals 
anomalously large finite-size effects at low frequencies while a frequency-moment analysis indicates a finite 
d.c. conductivity. This phenomenon also appears in a prototype integrable quantum system of impenetrable 
particles, representing a strong-coupling limit of both models. In the thermodynamic limit, the two results could 
converge to a finite d.c. conductivity rather than an ideal conductor or insulator scenario. 

PACS numbers: 71.27.+a, 71.10.Pm, 72.10.-d 



I. INTRODUCTION 

Transport of strongly interacting fermions in one- 
dimensional (ID) systems have been so far the subject of 
numerous theoretical as well as some experimental studies'. 
While the ground-state and low-temperature properties, fol- 
lowing the Luttinger-liquid universality, are well understood, 
the transport properties still lack some fundamental under- 
standing regarding the role of fermion correlations. It has 
become evident in recent years, that with respect to transport 
(in contrast to static quantities) integrable many-fermion mod- 
els behave very differently from nonintegrable onesifc^. Some 
basic ID fermion models are integrable, as the t-V model 
(equivalent to the anisotropic Heisenberg spin model) and the 
Hubbard model, and reveal in the metallic regime dissipation- 
less transport at finite temperature T > 0, well founded due to 
the relation to conserved quantities''*'^. The transport in the 
'insulating' regime of integrable models, however, has been 
controvertial and is the issue of this paper 

Let us concentrate on the dynamical conductivity in ID sys- 
tem 

a{uj) ~ 2'kD5{lj) -[^ Gregi'^), 
areg{^ > 0) = , Re / dte*'^* (j (0)) , (1) 
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where j is the (total) particle current operator, /3 — 1/T 
and L is the number of sites in the chain (we set everywhere 
fcs = S = eo = 1 as well as lattice spacing ap = 1). At finite 
T the charge stiffness (refeiTed to also as Drude weight) D{T) 
measures the dissipationless component in the response, while 
crreg(w) is the 'regular' part. The requirement that the ground 
state is insulating^ is £>o = D{T = Q) =0. In the insulat- 
ing regime there are still several alternative scenarios for the 
behavior at finite temperatures. The system can at T > be- 
have as: a) an 'ideal conductor' with D{T) > 0, b) a 'normal 
resistor' with D{T) = but cro = ^ 0) > 0, and c) an 
'ideal insulator' with D{T) =0 and cto = 0. 

A well-known T = insulator is the t-V model at half fill- 
ing and V/t > 2. The model is equivalent in this regime to an 



easy-axis anisotropic XXZ Heisenberg model with A > 1. It 
has been shown by one of the present authors^ that D{T > 0) 
is finite for V/t < 2 but decreasing towards D{T > 0) = 
at V/t — 2. This gives a strong indication that D{T) = in 
the whole regime V/t > 2, although there are also alterna- 
tive interpretations'^'^. The present authors speculated in this 
case on a possible realisation of an 'ideal insulator' ' where 
also a-o{T > 0) = 0. The argument is based on the observa- 
tion that at least in the V/t oo limit the soliton-antisoliton 
mapping can be applied, where the eigenstates cannot carry 
any current. However, the issue proved to be more involved. 
Note that the transport of gapped spin systems described by 
the quantum nonlinear sigma model, when treated by a semi- 
classical approach'*' (mapping to a model of classical impen- 
etrable particles) indicates a 'normal conductor' with a finite 
diffusion constant and cto(T > 0) > 0. On the other hand, 
a Bethe ansatz approachii concludes to a finite Drude weight 
D{T > 0) > and thus ballistic transport. It should be re- 
minded that the V/t — 2 case, corresponding to the most stud- 
ied isotropic Heisenberg model, is marginal situation, with the 
long-standing open question whether the diffusion constant 
(studied mostly at T ^ oo) in this model is finite^ '-. Another 
prominent T = insulator is the Hubbard model at half fill- 
ing. Here even the question of D{T > 0) is controvertial. On 
the basis of Bethe ansatz results''' and Quantum Monte Carlo 
simulations'** it is claimed that D{T > 0) > 0, i.e. an 'ideal 
conductor' situation. More recent analytical considerations- 
seem to favor D{T) — 0. 

The aim of this paper is to present numerical evidence that 
the dynamical conductivity fT(aj) in the insulating regime of 
several integrable ID models is indeed very anomalous. We 
consider in this context three ID models: the t-V model, the 
Hubbard model and a related model of impenetrable particles. 
First, finite-size scaling of results for all mentioned models 
indicates that indeed D{T > 0) ~ (whereby the evidence is 
somewhat less conclusive for the Hubbard model). Moreover, 
we show that on the one hand small-system results reveal large 
pseudogap features in cr(ti; 0) and large finite size effects 
extending to high frequencies; on the other hand, after the 
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finite-size scaling in the thermodynamic limit is performed, 
the results could be consistent with a 'normal' and featureless 
cr(w) found by a frequency-moment analysis. In this respect 
the behavior is very different from the one in nonintegrable 
quantum many -body models where even in small-size systems 
a 'normal' diffusive behavior is very evident^i^. 

The paper is organized as follows. In Sec. II we present two 
alternative numerical methods used to analyse the dynamical 
conductivity <t{llj): the microcanonical Lanczos method and 
the method of frequency moments. In Sec. Ill results for three 
different ID models in the insulating regime are presented and 
discussed: the t-V model at half-filling, the Hubbard model at 
half-filling and the model of impenetrable particles. 

II. NUMERICAL METHODS 

Microscopic models considered in this paper are ID tight- 
binding models with the hopping only among nearest neigh- 
bors. We investigate within these models the dynamical 
charge conductivity a{u!) (in the case of impenetrable parti- 
cles the related spin conductivity <Ts{uj)) at T ^ cx3 with an 
emphasis on the low cu ^ behavior. The first approach we 
apply is the full exact diagonalization (ED) of the Hamilto- 
nian on a lattice with L sites and periodic boundary conditions 
(p.b.c.) taking into account the number of fermions N and the 
wavevector q as good quantum numbers. E.g., this allows for 
an exact solution of (t{lu) up to i = 20 for the t-V model. 
Larger systems can be studied using the Lanczos method of 
ED. 

Particularly appropriate at large enough T is the micro- 
canonical Lanczos method (MCLM)'^. The MCLM uses the 
idea that dynamical autocorrelations (in a large enough sys- 
tem) can be evaluated with respect to a single wavefunction 
\'^) provided that the energy deviation 

(5e= A)2|^))i/2 (2) 

is small enough. Clearly, A determines here the temperature 
T for which j^*) is a relevant representative. Such \'^) can 
be generated via a first Lanczos procedure using instead of H 
a modified projection operator P — {H — Xf' , performing 
Ml Lanczos steps to get the ground state of P. The dynam- 
ical correlations are then calculated using the standard Lanc- 
zos procedure for dynamical autocorrelation functions, where 
the modified = il^*) is the starting wavefunction for the 
second Lanczos iteration with steps generating the con- 
tinued fraction representation of aiui). The main advantage 
of the MCLM is that it can reach systems equivalent in size to 
the usual ground-state calculations using the Lanczos method. 
For details we refer to Ref.[17]; e.g., the largest available size 
for the t-V model is thus L — 28. 

Besides the g[uj) spectra it is instructive to also show the 
normalized integrated intensity In tight-binding models 
with n.n. hopping the (optical) sum rule for oiuj) is given 
by {—T)/2L where T is the kinetic term in the Hamiltonian. 
Hence /(w) can be expressed as 

or 

I{uj) ^D* + j-^ j^^ doj'a{oj'), (3) 



which is monotonously increasing function with the limiting 
value I{u) oo) — 1 and well defined even for small sys- 
tems. Here, D* = 2LD/{-T). It should be noted that in a 
full ED calculation the Drude part D appears strictly at a; = 0, 
Eq. Q, while in the MCLM it spreads into a window Se ^ t 
governed by the number of Lanczos steps Mi . Still, choosing 
large enough Mi ^ 1000, Se becomes very small, hence we 
get well resolvable Drude contribution. Typically we use in 
the calculations presented here Mi = 1000, M2 = 5000. In 
order to get smooth spectra especially for small system sizes, 
we additionally performed an averaging over Nx different A 
with respect to the normal Gaussian distribution. Typically, 
we used N\ ^ 20 for smallest L and N\ ^ 1 for largest L 
presented in figures below. 

As will be evident from results furtheron (t{uj) exhibits 
huge finite-size effects. The latter are clearly a consequence 
of the integrability since nonintegrable models do not ex- 
hibit such features. In order to avoid such finite-size phe- 
nomena we also perform an alternative analysis using the 
method of frequency moments (FM). It is well known that at 
T — 00 one can calculate for (j{uj) exact frequency moments 
m2k = T^iJ-2k/T as 

M2fe = Tr([i/, [i/, • • • , [H,j]] ■ ■ •]i)/Tr(l). (4) 

Moments correspond here to an infinite system L 00 
and could be evaluated at fixed fermion concentration n — 
N/ L using the linked cluster expansion and the diagrammatic 
representation'** '^. Only clusters containing up to A: + 1 par- 
ticles can contribute to /Lt2fe in an infinite system for models 
with n.n. connections only. However, an analytic calculation 
of moments for larger k becomes very tedious. Hence we use 
the fact that exact moments for an infinite system can be ob- 
tained also via the ED results for small-system (with p.b.c. )2£ 
provided that the system size L is large enough. I.e., 

1 ^ 

N = 771,1 

where |vE'Arm) refer to eigenstates for N fermions. In Eq. (|5} 
^ = 'ENNst{N)f^ and fugacity / = exp(/i/T) can be 
related to the density 

''^^Y.^^^^(^)f''> (6) 

N 

where Ngt is the number of states for given N. Let us illustrate 
the feasibility of FM method for the ID t-V model. Perform- 
ing full ED for all fillings < < L on a ring we get exactly 
FM uptofc = L/2 — 1 whereby even higher k > L/2 — 1 
moments are quite accurate. Using full ED for i = 20 we 
thus reach for the t-V model exactly up to ^ig. 

Next step is to reconstruct spectra (t{lli) from ^2k with 
fc = 0, K. There are various strategies how to get the spectra 
most representative for K ^ 00, expecting a smooth function 
a{Lo). We follow here the procedure proposed by Nickel^-. 
First, a nonlinear transformation u — z + ('^ / z is performed 
where C, is chosen as the largest eigenvalue in a truncated con- 
tinued fraction representation of fT(aj). For (j{u!) then a Pade 
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approximant \Ki/K2\ is found in terms of functions of the 
novel variable z. 

It should, however, be noted that FM are less sensitive to 
the low-w regime, so a possible price to pay is an uncertainty 
in the low frequency results. In this respect MCLM and FM 
results really yield an alternative view of low-w dynamics. 



III. RESULTS 
A. t-V model 

Let us first analyse the \Ti t-V model for interacting spin- 
less fermions, 

H = -tY^{c\^^c^ + h.c.) + V^nini+i, (7) 

i i 

with the repulsion V between fermions on n.n. sites and the 
corresponding current operator 



ED: L=16 
ED: L=20 
MCLM: L=24 
MCLM: L=28 
FM: [4/5] 
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-t'^iicj^j^c. + h.c). 



(8) 




At half-filling, i.e., at the fermion density n = 1/2, the ground 
state is metallic for V/t < 2 and insulating for V/t > 2. 
Note that by introducing a fictitious magnetic flux via the sub- 
stitution t <e"^ the model turns into the anisotropic XXZ 
Heisenberg model for (j) — tt/L and even number of fermions. 

In the following we present only results in the limit T 
oo. From Eq. Q it follows that a scales in this limit as (3, 
hence we present in Fig. la a{uj) / (3, calculated for even num- 
ber of fermions for V/t — A and various sizes L = 16 — 28. 
Results for D/ (3 are plotted in the inset of Fig. lb and reveal 
an exponential decrease with L, which is at the same time a 
challenging test for the feasibility and the sensitivity of the 
MCLM at lai-ger sizes L. 

From results in Figs. la,b several observations follow: a) 
the dissipationless component D becomes negligible at large 
L and the extrapolated value for L = cx) is consistent with 
Z) = ^, b) there is a pseudogap at low oj followed by a pro- 
nounced peak at a; = and damped oscillatory features at 
bj > LUp, almost up to the bandwidth ~ At, c) the peak and ac- 
companied oscillations move downward with the system size 
approximately as ujp oc 1/L, d) the pseudogap in a{uj 0) 
is compensated by the peak intensity as evident from the inte- 
grated I{uj) which is essentially independent of L for ui > ujp, 
e) I{uj < t) could approach I{uj) ^ ctolu for large L, indica- 
tive of a 'normal' d.c. conductivity (Tq in the thermodynamic 
limit. 

When applying the FM method to the t-V model we get 

r! = (l-n)-^, f = n/{l-n). (9) 

Using full ED for L = 20 we reach exactly up to nig. In 
Figs, la.b we diplay results for (t{uj) obtained via the FM us- 
ing K = 9 and the corresponding [4/5] Pade approximant. 
The FM method proves to be very stable in particular with re- 
spect to the most interesting and sensitive value ctq. Namely 



Figure 1: a) Conductivity (7(aj)//3 and b) integrated normalized I{w) 
at T — > cx) witliin the ID t-V model witli V/t — 4, obtained using 
the ED and the MCLM for systems with length L, and the frequency 
moment expansion. The inset shows ln(Z)//3) vs. L, where the line 
is guide to the eye. 

the latter varies only slightly between, e.g., [3/3] and [4/5] 
Pade approximant. Results confirm the overall agreement of 
MCLM and FM-method spectra apart from evident finite-size 
phenomena at < ujp. It should be, however, mentioned that 
there are still some nonessential differences between I{uj) re- 
sults even at higher uj > ujp since the MCLM results are for 
fixed fermion number N = L/2 whereas the FM corresponds 
to a grandcanonical averaging over all N so that even lowest 
moments differ slightly. The general conclusion of the FM ap- 
proach is that it does not show any sign of pseudogap features 
and thus favors quite featureless a{oj) with finite ctq- Essen- 
tially the same results are reproduced for (t{uj) analysing FM 
using the maximum-entropy method-^. 



B. Hubbard model 

Next let us consider the ID Hubbard model 

H = -t^(e''^cJ^.i,Q, +h.c.) + 

i.s i 

j = -i^(ze^^4i^,c,,+h.c.), (10) 

i.s 

where we take into account a possible fictitious flux 0. We 
study the model at half filling n = N/L ~ 1 where the ground 
state is insulating, i.e. Dq = 0, for all L/ > 0. In the limit 
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L ^ oo the behavior should not depend on (p. Nevertheless 
in small systems low-o; features, in particular D{(p), depend 
on (j>. We present here calculations within the Hubbard model 
using the ED and the MCLM at = tt/ (2L) since in this case 
D{(/), T) is at maximum. Relative to the t-V model, smaller 
sizes are reachable for the Hubbard model at n = 1, i.e., we 
investigate L ~ 10 performing full ED, while with the MCLM 
systems up to i = 16 can be studied. 




Figure 2: a) cr(ci;) / f3 and b) I{w) within the ID Hubbard model with 
U/t = 4, obtained via the ED, the MCLM (finite L) and the FM 
method. The inset shows D/P scaled vs. 1/L, whereby the Une is 
guide to the eye. 

Results for the intermediate case U/t = 4 and again T 
oo are shown in Figs. 2a,b. We note that several features are 
similar to results for the t-V model: a) D decreases with L, 
b) a pseudogap appears for oj < ujp, c) large finite size effects 
extend up to frequencies of the order of the bandwidth, d) the 
pseudogap scale appears to close with the increasing system 
size. 

However, the dependence of D{L) is not exponential, but 
the scaling appears to follow D oc 1/L (see the inset of 
Fig. 2b). Although with less certainty than within the t-V 
model we could again support the limiting value D{T) — 0. 
Also, I{lu) tends with increasing L Xo I ^ (Tquj for oj < t, 
here approaching from higher values in contrast to Fig. lb. 
In spite of differences to the t-V model, results scaled to the 
thermodynamic Umit could be consistent with a smooth (t{lu) 
and a finite o-q- 

We also perform the FM analysis, using exact ED results 
for systems with up to i = 10 and < < 2L. Here, we 



The analysis is accurate up to /iio and corresponding [3/2] 
Pade approximants. This is barely enough to reproduce gross 
features of limiting a{uj), nevertheless results are in agree- 
ment with previous conclusions for the t-V model. 



C. Impenetrable particles 

The above results indicate that integrable models in the 'in- 
sulating' regime share similar features in the dynamical con- 
ductivity a{uj). It has already been proposed^ that it is helpful 
to consider the large interaction limit, i.e., V ^ t and U ^ t, 
where the dynamics of both models is simplified but remains 
highly nontrivial. For a half-filled band in this limit we are 
dealing with an excitation spectrum composed of split sub- 
spaces with a fixed number Ng of oppositely charged "soliton- 
antisoliton" (ss) pairs. In such a limit, the solitons/antisolitons 
- doubly occupied/empty sites in the Hubbard model, occu- 
pied/empty n.n. sites in the t-V model - behave as impen- 
etrable quantum particles, since their crossing would require 
virtual processes with AE = U, V. 

The simplest prototype model which incorporates the same 
physics - that of a system with two species of impenetrable 
particles - is the ID t-model. 



H = -i^(5|+i_^Ci^ +h.c.), 



(12) 



where projected fermion operators take into account that the 
double occupation of sites is forbidden; the two species of 
particles are represented by the up/down spin fermions. Thus 
we consider within the i-model the spin current 



i+l,s 



h.c), 



(13) 



2L 



/ = n/(2-n). 



(11) 



and the corresponding spin diffusivity as{oj). 

The only relevant parameter within the t-model is the elec- 
tron density n ~ n-^ + ni, where < n < 1 and of interest 
is the paramagnetic case = ni- The model ( I12t is also 
exactly solvable. Moreover, the electron current j commutes 
with H, while the spin current jg does not. It is plausible 
that in an unpolarized ring, — N^, exact eigenstates do 
not carry any spin current, i.e., {'^n\js\'^n) = 0, and hence 
D{T) = 0. This becomes clear by introducing the ficti- 
tious flux by t — > ie*"^. Particles cannot cross, so all eigen- 
ergies e„ are independent of (p. Since D{T) can be related^ 
to d'^en/dcfP this leads to D{T) = 0. Stifl, this does not pre- 
clude (Ts(w > 0) > 0, since (^'nljsl^'m) ^ in general. 

We study o's(w) within the t-model again using the same 
methods. With the full ED we reach L = 12 while with the 
MCLM up to L = 20 sites. For the presentation we choose 
the quarter-filled case, n — 1/2, where most systems are 
available, L = 8, 12, 16, 20. Results are shown in Figs. 3a,b. 
As expected, finite-size features are very similar to those in 
Figs. 1,2, apart from D = Q. The pseudogap is pronounced 
even more clearly, with the peak frequency LOp cx 1/L. Par- 
ticularly powerful for this model is the FM method. Namely, 
from the full ED results we can evaluate exactly moments up 
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to /i2o- Since there is a single characteristic scale t, the struc- 
ture of as{u!) is simpler and better reproducible via the FM 
method. Results corresponding to [5/5] Pade approximant are 
presented in Figs. 3a,b and again indicate on a 'normal' diffu- 
sivity in the thermodynamic limit. 




Figure 3: a) spin diffusivity as{uj)/f3 and b) the integrated 
within the ID i-model at n = 1 /2, obtained via the ED, the MCLM 
and the FM method. 



IV. CONCLUSIONS 

Let us summarize and comment obtained results. We have 
shown that all considered ID integrable models of interacting 
fermions share several common features: 



a) The charge stiffness is either D = (i-model) or appears 
to scale to zero, whereby the evidence is stronger for the t-V 
model. Whereas the vanishing D is easy to understand for im- 
penetrable particles, it is a highly nontrivial statement for the 
other two model s^i^i'^i'^i'^ . The observations can be rational- 
ized in a way that the t-V and Hubbard model at half-filling in 
the thermodynamic limit L ^ oo remain to behave as in the 
limit V,U ^ oo where solitons and antisolitons cannot cross. 

b) The pseudogap is pronounced for finite-size systems 
whereby the finite-size peak scales as uip oc 1/L. 

c) The extrapolation to the thermodynamic limit could be 
compatible with a rather featureless and regular (t{uj ^ 0) 
and thus finite cto- With respect to the last two points the 
ED (including MCLM) and FM methods are complementary. 
Whereas the FM method (valid for an L ^ oo system) cannot 
detect finite-size effects and appears to converge to a feature- 
less <t{llj), the ED methods are evidently sensitive to the effect 
of p.b.c. at finite L. 

A fundamental question raised by these observations is, 
whether the large finite size effects observed at low frequen- 
cies are reflected to the dynamics of bulk systems and in 
particular, which features of the conductivity (e.g. q,u! - 
dependence) might be singular 

We restricted our results to T ^ cxd. The latter is clearly 
most convenient for the FM method. Nevertheless, from the 
MCLM results considered at finite but high T there appears 
no evidence for any qualitative change on behavior We also 
presented results for a single parameter for the t-V and Hub- 
bard model, and one filling for the t model, respectively. One 
cannot expect any essential difference for other values within 
the 'insulating' regime, although numerical evidence becomes 
poorer, e.g., on approaching V ^ 2 t within the t-V model. 
Clearly, the most challenging case is V = 2 t, corresponding 
to the isotropic Heisenberg model. Our results reveal an in- 
crease of (To on approaching V ~ 2 t. Still we are not able to 
settle the well-known dilemma^' whether cro remains finite 
or diverges in this marginal case. 
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istry of Education, Science and Sports, under grant PI -0044. 
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